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^ ■ Abstract 
O . 

. We determine the mass dependent renormalization as well as improvement coefficients for the heavy- 
light vector and axial-vector currents consisting of the relativistic heavy and the domain-wall light 
■ quarks through the standard matching procedure. The calculation is carried out perturbatively at 

(N : 

the one loop level to remove the systematic error of 0{as{amQ)^ap) as well as 0{as{amQ)'^) (n >0), 
where p is a typical momentum scale in the heavy-light system. We point out that renormalization 
' and improvement coefficients of the heavy-light vector current agree with those of the axial-vector 



o 
o 



current, thanks to the exact chiral symmetry for the light quark. The results obtained with three 
different gauge actions, plaquette, Iwasaki and DBW2, are presented as a function of heavy quark 
\ mass and domain- wall height. 
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I. INTRODUCTION 



The study of heavy-light hadrons including a single heavy quark {b or c) has played a central 
role in the particle physics. In particular, much attentions have been paid to estimating weak 
decay matrix elements of these hadrons to determine the Cabbibo-Kobayashi-Maskawa matrix 
elements. Lattice QCD simulation offers a promising method to calculate the decay amplitudes 
from the first principles, and a number of efforts using lattice QCD have been already made 
and contributed to the heavy-light physics 

However, the lattice calculation of the heavy-light system is not yet precise enough. Since 
the currently accessible inverse lattice spacing is at most 1/a ~3, 4 GeV, so that arrive 
never holds, the lattice calculations suffer from sizable 0((amfec)") {n >1) errors. Another 
difficulty, though this is a rather general problem in Lattice QCD simulations of hadrons, arises 
from chiral extrapolations of physical quantities to their physical values (~ as numeri- 

cally investigated in Refs. 3, 0, Q . Because of limited computational resources, we are forced 
to simulate light quarks much heavier than the physical light quark masses, and, then, long 
extrapolations are required. In addition, non-analytic behaviors of physical quantities around 
the massless point which chiral perturbation theory predicts make reliable extrapolations even 
more difficult. 

Onepossible way to overcome the above obstacles is to combine the relativistic heavy 
quark 6| and the domain- wall light quark to simulate heavy-light quark systems. 

The relativistic heavy quark [o] is obtained by pushing the on-shell improvement pro gram [lo| 
to the Wilson type massive fermions, and its action is equal to the Fermilab action ll| with 
a special choice of parameters. The relativistic heavy quark action has four parameters to be 
tuned, and the perturbative determination of these parameters as well as the wave function 
renormalization has been completed to the one-loop level so far 12] . As for the problem of 
chiral extrapolations in terms of light quark masses, we employ the domain-wall formalism for 
light quarks, which has already been used successfully in large scale simulations. It has been 
demonstrated that this formalism allows to simulate lighter quarks than the standard Wilson 
type quark action. 

In this paper, we determine the renormalization coefficients for the heavy-light vector and 
axial-vector currents consisting of the relativistic heavy and the domain-wall light quarks 
through the standard one-loop matching. The improvement coefficients associated with the 
dimension four operators are evaluated to the one-loop as well. Since the obtained coefficients 
depend on both heavy quark mass tuq and domain- wall height M5, we work out the interpola- 
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tion formula to express the results as simple as possible. In this work we take three different 
gauge actions, the plaquette, the Iwasaki and DBW2 14, 3|; as examples. The Feynman 
rules and useful formula for the gauge actions are described in Ref. in rather general man- 
ner, so we will omit that part of explanations from this paper. With the coefficients obtained 
here and the lattice calculations of all relevant operators, we can remove the systematic error 
of 0{as{amQ)"'ap) as well as OlaglamQ)") {n >0) from physical amplitudes. As for the cal- 
culational method, we completely follows that in Ref. [16|, where the similar calculations are 
performed using the clover action for the light quark. 

This paper is organized as follows. In Sec. |nj we briefly explain the quark actions used 
in the following sections. In Sec lIIIl the definitions of the renormalization and improvement 
coefficients which we determine in this paper are given. Then, we give a short notice on the 
benefit of chiral symmetry in Sec JIVl The calculational method and the numerical results 
are given in Sees. In Sec. IVH we consider two ways to improve the convergence of the 
perturbative series. Concluding remarks are stated in Sec. lVIII In the process of this calculation, 
we found that the on-shell improvement of the massive domain-wall fermions is necessary. The 
improvement procedure in this case is briefiy explained in appendix El 



II. THE QUARK ACTIONS 



A. The relativistic heavy quark 



In this work, the heavy quark is described by the relativistic heavy quark action 



0, 



x) 



Sq = Y, moQ{x)Q{x) + Q{x)'^qDqQ{x) + vqY^ Q{x)'jiDiQ{ 

X L i 

--fQ{x)DlQ{x) -'-fY.Qi^)DlQ{x) 

--^CE'^Q{x)(JoiFoiQ{x) - ^-^CB^Q{x)aijFijQ{x) 

The parameters appearing in the action are determined by requiring n-point amputated on- 
shell Green's function on the lattice to be equal to the continuum one up to O(a^). In order to 
obtain the improvement coefficients which are valid at an arbitrary quark mass, no expansion 
in (amo) nor l/(amo) is allowed through the calculation. By matching the two-point function 
(propagator) to the continuum one at tree level, the pole mass, wave function, z/g and are 
found to be 
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Matching the amputated three-point function (quark-gluon form factor) leads to 
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Since Tt is not constrained from the improvement condition, in the following we set rj=l for 
simplicity. Once the parameters have been tuned to their proper values given in eqs. (0)^©; 
we can remove the systematic errors of ©((amg)") and 0{{amQ)^ap) {n >0) from the action. 
Extending the above procedure to the one-loop level, O(ajamo)"') and 0{as{amQ)"'ap) 



12], and the numerical 



(n >0) errors are removed. This calculation was carried out in Ref. 
results depending on the quark mass are given in terms of an interpolation formula. Here we 
only quote the result of the wave function renormalization for the later use. Defining the wave 
function renormalization by 
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AQ{mQ ) is, to good precision, given by 



Aq(0) + 
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121. 



where the values of Aq(0), Oj and 6j (i = 1, . . . , 4) are listed in Table II of Ref. 

It should be noted that in the massless limit all t'g'', r^^\ c^^ and c^^ are reduced to unity (for 
rt=l), and the action is reduced to the 0{a) improved SW quark action [17]. More importantly, 
the space-time rotational symmetry is restored in this limit to all orders of a and a^- 



B. The domain- wall quark action 

As for light quarks, we take the domain-wall fermion 

BBS 

. More precisely, we employ 

the following version of the domain-wall action in this work, 

5'dw = XlZlV^''(^)^Dw(a;, s;i/,t)^t(l/), (11) 

x,y s,t 



DD^^{x,s;y,t) = Dy^{x,y) Ss,t + 5x,yDr,{s,t), (12) 
Dw{x, = E 7m^m - f E - ^5 (13) 

D5{s,t) = Ss,t - PlSs+l^t - PBSs-l,t + mf[PL6s,NsKt + PRSs,lSN5,t] , (14) 

where Pr/l = (1 i 75)/2, and s and t label the coordinate of the fifth dimension running 1 to 
N^. We set A^s ^ oo throughout this paper. The physical quark field q{x) is, then, constructed 

by 

q{x) = Piipiix) + Pni^Nsix), (15) 
q{x) = Mx)Pn + ^/jN,{x)PL. (16) 

As long as massless quarks are concerned, this original action is enough to proceed. However, as 
explained in Ref. [l6[ and later as well, in a part of calculations we should keep the light quark 
mass small but finite. Once quark masses become finite, the infrared divergences which do not 
exist in the continuum show up in a certain class of loop integrations due to the mismatch of 
dispersion relation and make perturbative matching to 0{asap) meaningless. To eliminate the 
mismatch, the original domain-wall action is required to be improved, and the improvement 
can be realized by simply applying the on-shell improvement program P, Q| to the massive 
domain- wall fermion. As a consequence, it is found that we need, at least, two modifications 
in the original domain-wall action. First we replace the first term of eq. ()13|) with 



and, then, add a new term, 

- Y Qix) E A' Qix) = 'Y ^s{x) E A' (Pl Ss,nAi + Pr 5sA,n,) ^*(x), (18) 

i i 

to the action eq. (fTT|) . Here we introduced two parameters Vq and Rg, which can be tuned by 
following the on-shell improvement program [^. As a result (see appendix^for the derivations), 
we obtain 

sinh(m^°)) ^ ^ 

= Tsf^. (19) 
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where m^°) is the tree- level pole mass depending on and M5, and its expression is given in 
eq. ()A8|) . As nif ^ 0, these parameters approach to 

^f-l + 0(mj), (21) 
^ ^ + 0(mj), (22) 

respectively, and, thus, the action is reduced to the original one. Analytical studies of the 
loop integrations reveal that the tree level matching of Ug and Rs suffices to get rid of the extra 
infrared divergences associated with massive light quarks, so that it is not necessary to introduce 
terms like the clover terms in this work though to remove the 0{{amf^)''^ap) errors from the 
action throughly we still need the clover terms whose tree-level coefficients are calculated in 
appendix 1X1 

For the wave function renormalization, the one for the massless domain-wall fermions is 
sufficient for the following calculation. The relation of the massless domain-wall quark field to 
the continuum one with MS in NDR is given by 
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Wo = 1- Ms, (24) 
Zgiwo) = 1 + 9' (-^ In(aV') + A,(w;o)) , (25) 
Z^iwo) = l + g^A^{wo), (26) 



and the numerical values of Ag{wo) and A^(wo) for various values of wq 
from Z2(DRED) and S^;, which are tabulated in Tables I, III and V of Ref. 
Iwasaki and DBW2 gauge actions respectively, through the relations. 



M5) are constructed 
for the plaquette. 



^.(^0) = ^;^2(NDR) = (^2(DRED) + 1) , (27) 



A.(^o) = T^S- (28) 



Cf 2wo 
167r2 1 — w'q 

The space-time rotational symmetry is recovered in the massless limit to all orders of a and a^. 



III. DEFINITION OF THE IMPROVEMENT COEFFICIENTS 

Taking into account existing symmetries such as parity, charge conjugation, etc., it is turned 
out that the renormalized continuum vector and axial-vector currents can be expressed in terms 
of 0(a) improved lattice operators as 



T^MS/^^ _ vMS-latt 



(l{xh^,Q{x) - g^c^^d^ {q{x)Q{x)} - g^Cy^d^ {q{x)Q{x)} 



-g'4^{d,q{x)HjM^) - 9'clq{x)j,j,{d,Q{x)} + 0{g')\ , (29) 
q{xhf,i5Q{x) - g'^cX^d^{q{x)-f5Q{x)} - g'^c^^d^{q{x)-f5Q{x)} 
+g'4Mqix)HwQix) - /cf^g(x)7^757i{^.Q(x)} + Oig^)] .(30) 



where d^=d^ ^d^, q{x) and Q{x) represent light and heavy quark fields defined on the lattice, 
respectively, and the coefficients Zy^^^^^^^ and Cy^f^^^ are the ones we determine in the following 
sections. While the right hand side can be written in a gauge invariant way by replacing partial 
derivatives with covariant derivatives, the above expression suffices for the following discussion. 
Although the general form given in eqs. ()29p and (jHUj) does not depend on details of the lattice 
actions, the coefficients Z^^^^^^^ and Cy^^^ do. Since the space-time rotational symmetry is, 
in general, violating on the lattice unless quarks are massless, the fourth and fifth terms on 
the right hand side are required to improve the violation. Therefore, although the equation of 
motion always allows us to set Cy^ = Cy^ Aq ~ ^ take this, it never means that the 

space-time rotational symmetry is recovered in the temporal component of the currents. 

IV. BENEFIT OF CHIRAL SYMMETRY 

Here let us mention an advantage of using lattice chiral fermions to describe the light quark in 
heavy-light system by tak,„g the he»y-Ught vector and a.ia.-veetor currents as examples. The 
similar discussion is made in Ref. |19j], where the renormalization of the static-light current is 
calculated using the overlap fermion for the light quark. The authors of Ref. also mentioned 
the similar point in the calculation with the NRQCD heavy and the naive light quarks. 

Now let us consider a system including Nf=2 massless light quarks g*(x) {i = 1,2), de- 
scribed by the lattice fermions with the exact chiral symmetry, and one heavy quark Q{x). An 
infinitesimal fiavor non-singlet axial transformation is defined by 

6Aq'{x) = iw^{x)^^j,q^{x), (31) 

6A(f{x) = q^{x)iw\x)^^^^^, (32) 

where i and j are the fiavor indexes and run 1 to A^/, and is the Pauli matrix. Hereafter we 
take an arbitrary transformation parameter w"'{x) as 

{5a 3 for X G R 
(33) 
otherwise 

where i? is a certain space-time region with smooth boundary DR. The exact chiral symmetry 
for the light quarks ensures the existence of the conserved axial- vector current. Using this fact. 



we can derive the following axial Ward-Takahashi identity on the lattice 

(Oh / da,{y) Af\y) O,^, O^xt |0) = (0| (^A^int) ^ext |0), (34) 

JdR 

where A^'^^(?/) is the third component of the lattice conserved light-light flavor non-singlet axial- 
vector current, and Omt and Oext are arbitrary operators localized in the interior and exterior 
of i?, respectively. Now suppose 

0,^, = V^''-\x)=(t{x)^Mx). (35) 

then, 

^^Oint = -iqK^)\{ry'lul,Q{x) = -^(r3)^Mf '(x). (36) 
is read. Substituting the above into eq. (jSH), one would reach 



(0| l ja,{y)A^f{y) V^'^^^ix) Oext |0) = -- (0| A];'^\x) Oext |0), (37) 



1 

for i = \. 

Now repeat the above discussion with the renormalized currents, and rewriting them in 
terms of the bare currents, one obtains 

^MS-iatt^O| /• i^^[y)Af{y) Vl^\x) axt |0) = -Z^'^"^ (0| ^xt |0), (38) 

Here it is important to note that the conserved current does not receive the renormalization. 
Comparing eq. (|H7j) to eq. (|HHj). the relation 



^•MS-latt _ ^■'^S-latt /oq\ 
'V^i All V ' 



is obtained. The operator Oi^t is arbitrary, so that repeating the above procedure with the 
right hand side of eq.f^^ and comparing the resultant expression with eq. (jHUI) . 

4, = -^A,^ (40) 

cy, = (41) 

< = cl, (42) 

4 = (43) 

are derived in addition to eq. ()39|1 . Thus it is proved that, whatever the heavy quark action is, 
if a light quark field on the lattice has chiral symmetry, the renormalization constants and the 
improvement coefficients for vector and axial-vector currents are related to each other. This 
statement can be extended to any pair of operators which belong to the same chiral multiplet. 
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V. ONE-LOOP DETERMINATION OF THE COEFFICIENTS 



A. Calculational method 

In this section, the calculational method is outlined by taking the vector current as an 
example. More complete explanations are found in Ref. [l6|, where the similar calculation is 
performed with the clover light quark. Let us first define the off-shell vertex function of the 
lattice vector current q{x)'y^Q{x) by 

A'^'\p,q,mQ,m,,M,) = Y.{gjKf''^'\p,q,mQ,m,,M,). (44) 

The tree level contribution A^°-'=7^j is immediately followed. The one-loop contribution 
A^^\p, q,mQ,mg, M^), obtained by calculating Fig. ^ can be written in the following general 
form, 

+4ljF^ + 4iJsFj^ + i,-ijF^ + iJsl^F^ + iisl.F^^ 
+ {V, + q,) [G'^ + i>G^2 + + mi] 

+ {p, - q,) + m + 4^ + + 0{a'), (45) 

where 

3 3 3 3 

^ = Y1 P«7«, 4 = Y1 ?"7a, = J^Pi'^i^ is = J2 ^»7i, (46) 

a=0 a=0 i=l i=l 

and the coefficients F/, Gj* and Hj are the scalar functions of p^, q"^, p ■ q, ttiq, niq and M5. 
Here it should be noted that in fi=0 case we can set Fj=0 for /=3, 5, 7, 8, 9 and 10 without 
loss of generality. 

Sandwiching eq. (jl^j) by the on-shell quark states uqip) and Ug{q), which satisfy ]I>uq{p) = 
imgUQ^p) and Uq{q)4 = inigUg^q), the matrix elements are reduced to 

«,(g)A«'''^"(p, q, mg, m„ M5)uq{p) 
= X';'%{q)^^UQip) + Rf'ug{q)^^isUQ{p) + S'^''ug{q) isl^Uqip) 

+ Y';'\p, + q,)u,{q)uQ{p) + Z'^'\p, - q,)u,{q)uQ{p) + 0{a^), (47) 

where we defined 

^latt ^ + + ^^^p^ _ mgrrigF,^, (48) 

yutt = G^,+ irngG^ + im^G^^ - mQm.Gl, (49) 
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(50) 
(51) 
(52) 



By repeating the above procedure in the continuum, we obtain the continuum counterparts, 

^cont = Q in continuum because of the presence of the 



^T'^ y™''* and Z™'^*. 



space-time rotational symmetry. For X^""^*, the ultraviolet divergence is subtracted by the MS 
scheme, and using the resulting expression we define X^""^*'^ as 



■V^cont 

" 167r2 



+ ln(47r) -7s) = + 



(53) 



where we factored out the ln(/i2) term in order to make the /i dependence clear. The renormal 
ization factor of the vector currents is then given by 
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where 



Xlatt \r 



^7m 



cont,R 



(56) 



The ln(a^/i^) term in the first line of eq. is canceled by those in eqs. and (j2Sl), so 
the renormalization factor Zy^~^^^^ is independent of the renormalization scale /i. Although we 
take the naive dimensional regularization (NDR) throughout this paper to regularize ultraviolet 
divergences, the conversion to the dimensional reduction (DRED) can be made by 

1 



for both vector and axial-vector currents. 

The coefficients Cy'^'^ defied in eq. are found to be given by 
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Since the space-time rotational symmetry is restored at mq 



rrin 



even on the lattice, 



R^k^^\mQ=mq=0 — >S'Jf**'|m,Q=mg=0 — 0, hcnCC, Cy^ \mQ=mq=0 — 0. 

±,H,L 



'fc \mQ=mq=0 



To calculate the coefficients A^^ and c 



we define the difference between the lattice 



one-loop vertex function and the continuum one by 

A«(p, q, mg, m„ M5) = A«''^**(p, q, mg, m„ M5) - A«''=°'^Hp, g, mg, m„ M5). (62) 

While each of A^^^"''^** and A^^^''^"^^ has infrared divergences, they completely cancel in the 
difference. After subtracting the ultraviolet divergences and yU dependent terms from A^^^'^°^^, 
A^^ and cp^^'^ are extracted by proper 7-matrix projections on A^^^^ as follows, 
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-Tr 
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dX 



(1) 



1 + 70) 



;i+7o) 
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1. 

4 
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-Tr 
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X^f}\l + 7o)7fc + 2zm,^^(l + 70)7^7^ 



dX 
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dX. 
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dp. 



dqk 



p=p ,q=q^ 
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-Tr 

4 



Xi^\l + 70) + 2tmg—^{l + 7o)7fc 



p=p*^q=q* 



A. 



-Tr 



Af^(l + 7o) 



p=p ,q=qi 



(63) 
(64) 

(65) 
(66) 
(67) 
(68) 

(69) 
(70) 



where i, fc=l,2,3 and the external momentum p is set to p* = {imq, 0) and g to g* = {irriq, 0) 
or 12 = i-img,6). 

All coefficients A^^ and Cy'^'^ are obtained from the proper linear combinations of eqs. ()63p - 
(fTOj). As seen from eqs. (jU7j) . and (fTOj) . however, Cy^ and Cy^ can not be extracted at mg'^°)=0 
directly. Therefore, we keep the light quark mass small but finite in the calculations of these 
coefficients, where the on-shell improvement for the massive domain-wall fermions described 
in Sec. IIIII and appendix El is required to remove the associated infrared divergences. We 
estimate these coefficients with two different values of light quark mass mq^°)=0.0001, 0.001, 
and extrapolated linearly to their massless limit. 
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The loop integrations on the lattice are performed by a mode sum for a periodic box of 
a size with L=48 after transforming the momentum variable fc^ to k'^ = — sin/c^. To 
search for the mq and M5 dependences, we take 18 or more values of mg, depending on M5, 
in 0.1 < mg'^^) < 10.0 and 11 values of M5 in 0.1 < M5 < 1.9. We choose three different 
gauge actions, the plaquette, Iwasaki and DBW2, in this work. To see the uncertainties, a 



large part of integrations are performed with BASES [21 1 as well. For almost all of rrLq and 
M5 that we investigated, all the coefficients using BASES show a nice agreement with those 
with the mode sum, while the visible discrepancies appear in A^^, Cy^ and Cy^ at M5=0.1 and 
1.9, and the discrepancies grow as mq increases. Repeating the corresponding integrals using 
the mode sum with several values of L, it is found that the results approach to the BASES 
results as L increases. From the above observations, we restrict the data within rrLq < 7.0 and 
0.1 < M5 < 1.9. After this restriction, the remaining largest discrepancy reduces to less than 
0.001, which can only affect physical amplitudes by 0.5%, at most. 

The continuum part is evaluated in the naive dimensional regularization (NDR) with the 
modified minimal subtraction scheme (MS). To make cancellations of infrared divergences 
explicit suitable integrands, which can be estimated analyticallv, are subtracted from the lattice 
integrands and added to the continuum ones as done in Ref. 



B. Numerical results 

Figures EE show the niq'^^^ dependence of A^^^, Cy^^, Cy^^ and Cy^ , respectively. In each 
figure the results with the plaquette, the Iwasaki and DBW2 gauge actions are laid from top 
to bottom. We pick up the results with M5=0.5, 1.1 and 1.7 to give some idea about the M5 
dependence, which is turned out to be smaller than the mq dependence in the most cases. The 
solid lines in each figure are obtained by fitting the numerical data to 



where C, A^^ {jn=l,. . .,4) and D„(n=l,. . .,4) are again parameterized by 

C = ^c«(M5r, (72) 

i=0 

= En«(M5)\ (73) 

i=0 

Dn = jZd^HM.y. (74) 



i=0 
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45 parameters are thus determined by the fit. Their numerical values are tabulated in Tabs. |1J- 
IVIIII for each improvement coefficient with three different gauge actions. The ffi lines reproduce 
the data points very well over the whole region of ttiq' and M5 as seen in the figures. 

In the limit of niQ^^^ = mg^^^=0, the presence of space-time rotational symmetry guarantees 
that Cy^^ = independently of the gauge action and M5. So we put the constraint of c*^*^=0 
(2=0,1,2,3,4) in the ffi of these coefficients. The space-time rotational symmetry also guarantees 
that A^^ = A^g and Cy^ = Cy^ at mg^^^ = mg(°) =0. These relations are not used in the ffis as 
constraints, but they are found to be satisfied within 1 % in any cases. 

We also calculated the renormalization and improvement coefficients for the axial-vector 
currents independently, and conffimed that the relations derived in eqs. (jH^ - ipSj) well hold 
within the expected numerical accuracy. 



VI. IMPROVING PERTURBATIVE SERIES 



A. mean field improvement 



Applying the mean field improvement |22] is expected to accelerate convergence of perturba- 
tive series. In this section we examine the mean field improvement using the results obtained in 



the previous section. Following Ref. 



12j |. the improvement of the wave function renormalization 



for the relativistic heavy quark is implemented by replacing eq. (jHI) with 



where 



u 



MC ryiO) / ~ (0) 



Matt' 



1 + 9' 



Ha /i ) + AQ(m^O + -Y^MF + g^^(o) ^^Q 



~ (0) 



(0) 



(0) 



sinh(m§^) + rt{cos\i{ni§) - I) + 3(rs(mg^) - 1) 



cosh (mg ) -|- Tt sinh(mQ 



^Cf sinh(mQ(o)) + ri(cosh(mQ(o)) - 1) + 3(r,(mQ(o)) - 1) 4 
I cosh [vtiQ ^^'1 ) -|- Tt smh [mQ^^'i) 



mQ^'>^-g''AmQ + Oig' 



(76) 



(T5) 



and is the mean plaquette value measured in a Monte Carlo simulation whose perturbative 
expression is 



.PT 



W ^ = 1- g^^TMF, 



(77) 



where Tmf= 1/8) 0.0525664 and 0.0191580 for the plaquette, the Iwasaki and DBW2 action, 
respectively. The improved wave function renormalization for the the domain-wall fermions is 
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given by [1 



u 



MC 



-|^ln(aV) + A,(M5) + ^TMF 



(1 - 



(7J 



where = 1 - = 1 - (M5 - 4(1 - m)). 

Now taking into account the above reorganizations of the wave function renormahzations, 
tuq^^^ and M5, the mean field improved renormahzation factor is given by 



7MS-latt,MF 



u 



MC 



\ (i-*o)^r(^o) 



MF 



1 - /At?: 



Tsr^"'^- 



2 dniQ^ 



(79) 
(80) 



B. quasi non-perturbative improvement 

Let Z^^~y^^^ be the overall renormahzation factor for the lattice heavy-heavy vector current 
Q'joQ- Similarly define Z^^y^^^^ for the light-light vector current q'joq. These factors can 
be determined non-perturbatively without any special devises. We can, then, reorganize the 

perturbative series eq. (jH^ by using one or both of these factors as done in Ref. [2^. For 

7MS-: 

'qq,V 

improvement, eq. ()54p is modified as 



example, if one uses the non-perturbative value of Z^^y ^'^^^'^^ together with the mean field 



^MS-latt,QNP.MP ^ zg^^^m 



~ (0)\ ^MS-latt.NP, 



'gqy 



Wo) 



1 - /A^^P+^^ 



^QNP+MF 



Av., (m. 



^A^L(^o) _ ^Tmf 



l^JnZgL 

2 SmgW ' 



where the perturbative expression 



(81) 
(82) 



^S-^^"M = l + /A^L(t.o) 



(83) 



is used, and Ay (wq) is given by 



A^^(«;o) 



-^(^y/A(DRED)-- 



^^^,-../.(NDR) . 16.^ 2 



B4) 



The numerical values of zy/yi(DRED) are available in Tables I, III and V of Ref. This 
way of the improvement has an advantage that the factor Zw{wq) is automatically included in 



^MS-latt,NP/ \ 



We plot Ay„, A^?/ and A^J^^+^^ in Fig. El where the results with M5=0.5, 1.1 and 1.7 for 
three gauge actions are shown as before. From the figure, it is seen that in the non-improved 
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case the size of the one-loop coefficient depends on the gauge action while once one introduced 
the quasi non-perturbative and mean field improvements together the dependence becomes 
less transparent. It is also interesting to see that the mean field improvement makes the mq 
dependence smaller. The similar findings are made for the improvements of Ay^ as well. 

VII. CONCLUSION AND DISCUSSION 

In this paper, we determined the renormalization coefficients for the heavy-light vector and 
axial- vector currents consisting of the relativistic heavy and the domain-wall light quarks. The 
improvement coefficients associated with the dimension four operators are also evaluated, so 
that the leading errors are reduced to O(a^) and 0{a?p'^). The results are presented as an 
interpolating formula with two arguments, the heavy quark mass mq and the domain-wall 
height M5. 

Using light quark actions with the exact chiral symmetry brings the useful relations between 
vector and axial- vector current renormalization/improvement coefficients. One of the relations 
^Ms-iatt _ ^Ms-iatt especially useful for the test of the soft pion relation f {ql,^^)= f b / U 
where f^{q^) is one of the form factors defined in B ^ nlu transition amplitude. In this test, 
the ratio of z^^~^^^^ / Z^^'^'"^^^ plays an important role, and the unknown two loop coefficient in 
the ratio has prevented the test from being more precise. So once one determined Cy'^'^ at the 
one loop level, the leading uncertainty in the test is reduced to O(a^ap). Furthermore the use 
of the lattice chiral fermions for the light quark is expected to make the chiral extrapolation 
easier than that of the Wilson type light quark. With the lattice chiral fermions, thus, one can 
expect the significant reduction of the systematic error in the test. 
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APPENDIX A: ON-SHELL IMPROVEMENT OF THE MASSIVE DOMAIN-WALL 
FERMIONS: TREE-LEVEL ANALYSIS 



First let us define the following quantities, 

^ = 7o sin(po) + Vq^li sin(pi 



P 



W{p) = l-M5 + 2^ 



sm 



2 



cosh (a) 



1 + W^{p) + p 



2\W{p)\ ' 
mf{p) = uif + 2 Rs^"'-"'^ 

i 

uj = 2-M5. 



„ , Pi 

sm I- 



(Al) 
(A2) 

(A3) 

(A4) 

(A5) 
(A6) 



With the modifications made in eqs. ()17|) and ()18|) . the physical quark propagator with mo- 
mentum p is given by 

-i^ + mf{p) ( 1 - W{p) e"") 



S,ip) 



(AT) 



( W{p) e° - 1 ) + m^(p) ( 1 - e-°) ' 
The pole mass is obtained by finding the position of pole in the complex po plane with Pi=0 
and given by 

-mfu'^ + ^{1 + ■mjy + m'j uj'^{lu'^ — 4) 



m(°) = In 



1 + m'j — 2711 fUj 
In the case oi rrif <^1, the pole mass reduces to 

= M5(2-M5)m/ + 0( 



m. 



(AJ 



(A9) 



Note that the correction starts from Olmj) relative to the leading one. This is consistent with 
the discussion of chiral symmetry made in Ref. [6|. Solving this inversely, the bare mass ruf 
can be written as a function of the pole mass as 



UJ (2 cosh(m^°)) — u) — sJuP' (2 cosh(mW) _ cu)2 - 4 (sinh(m(o)))^ 



2 sinh(m 



{0)> 



Defining the wave function renormalization of the DW quark field Zq by 



its tree level value. 



7(0) 

^(j.latt 



mj (1 + w?f) cosh(m^°)) — 2m^ sinh(m^°)) 



[1 — mj) sinh(mg 
16 



(0)> 



(AlO) 



(All) 



(A12) 



is obtained by calculating residue of propagator at po = im^\ In the limit of m/ <^1, a 
well-known factor 



7(0) 



rrii 



1 



40) M,{2-M,) 



+ 0{m}) 



(A13) 



is recovered. 



(0) and are determined by imposing the correct dispersion relation on the quark prop- 



agator. From the denominator, it is found that uj^^^ and R^^^ must satisfy 



(z/W)2 [l + m}-2mf sinh(mf) ) + i?f (smh(mf))^ ^ ^ _ 



sinh(m^*^)^ 



(0) 



1 + TTij ) cosh(m^°^) — 2 mj sinh(m^°'') +2m/ (^cj — cosh(m^°'*)^ (A14) 



while from the numerator we obtain eq. (fT9|l . Substituting eq. (fT9|l to eq. ()A14|1 . eq. (j20|l is 
derived. 

We can further extend the improvement by adding the clover terms to the action, which 
removes the 0{{amq)"'ap) errors thoroughly, though this is not necessary in the present work of 
the main text. Although there are several possible ways to include the clover terms, we choose 
to add the following form. 



-q{x) 



ig ^ ^ ^ ^ ig 



q{x) 



to the action eq. (|TT| . To tune the new parameters CE,q and Cb^i it is enough to consider a 
quark form factor scattering off a single gluon 1^. The tree level form factor -F^(pi,P2) is, 
then, given by 



FS{pi,P2)= -ig {Z, 



(0) N-1 

g,latt/ 



nif cosh(m^°'') ( 1 + m?f — 2m f sinh( 



(0)^ 



sinh(mg'''') (1 — m'j) 
UqiPihoT''Uq{p2) 



2mj (^co" — cosh(m^°))) C 



1 — m'j 



+ 



(0) 
E,q 



C 



(0) 
E,q 



Uq 



sinh(m^ 



(0)^ 



sinh(mf ) M,(pi)T%(p2) 



(A16) 



for /i = 0, and 



(0) N-l 



latt/l 



X 



v^^ m/ ( 1 -|- — 2m/ sinh(m^°''^ 



sinh(m, 
17 



) (1 - 4) 



+ 



C 



(0) 



V, 



(0) 



sinh(m 



(0)^ 



mj (^uj — cosh(m^ 



+ i{Pu+P2i)\ — 



(0)- 

q • 



sinh(m. 



{0)^ 



XUg{pi)T''Ug{p2 



(A17) 



for n = i. Now using eq. (fT^. (^111) and ()A12|) . it is turned out that setting 
2mj fcj — cosh(m^°^' 



rriq^ (1 — m?j) 

9N / w sinh(m^°^)\ ^ 9 . , / / sinh(m^°)) 

m/(l + m^) cosh(m^°)) ^q.^ ~ ^ m^sinh(m^°)) u; - ^ 



(0) 



m 



m 



W 



m^g^ sinh(mg°'') (1 — mV 



(A18) 



reproduce the correct form factor, 



K{pi,P2) = -^9 M,(pi)7^r%(p2) + 0(0X2)- 



(A20) 



Note that and C^'g vanish in the massless hmit, which means that the action does not 
have 0{a) errors in this hmit. This is consistent with the presence of chiral symmetry in the 
domain- wall fermion. 
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i 





1 


2 


3 


4 








plaquette 






c(») 


4.50349e-02 


6.08797e-04 


-6.18501e-04 


2.789996-05 


9.798886-05 


(») 
n\' 


2.19429e-03 


-1.16680e-02 


-4.02314e-03 


3.756246-02 


-1.627736-02 


(i) 

n\' 


3.20289e-03 


3.14583e-02 


-1.83939e-01 


1.226586-01 


-3.124926-02 


n\' 


1.06077e-02 


8.31176e-02 


2.83741e-02 


-1.18936e-02 


-7.778386-03 


n\' 


3.89916e-04 


6.63305e-02 


-2.93231e-02 


2.145126-02 


-4.351426-03 


df 


4.88310e+00 


4.68459e+01 


-1.26499e+01 


8.524406+00 


-3.881396+00 


df 


4.33555e+00 


4.37762e+01 


9.61199e+00 


-1.430736+01 


2.714296+00 


df 


1.53503e+00 


3.98540C+01 


-6.97580e+00 


3.70893e-01 


-2.30419C-01 


df 


4.13266e-02 


7.86025e+00 


-4.21872e+00 


2.144466+00 


-4.476146-01 








Iwasaki 






c« 


4.48266e-02 


3.98377e-04 


-4.08681e-04 


8.388756-05 


3.255046-05 




2.01108e-03 


-1.03944e-02 


6.00284e-03 


8.165156-03 


-4.661346-03 


(i) 


5.10189e-03 


4.09113e-02 


-1.10778e-01 


6.816776-02 


-1.67361e-02 




5.50012e-03 


6.34936e-02 


6.54794e-02 


-4.654066-02 


3.991996-03 


(*) 
n\' 


-3.75123e-04 


2.81289e-02 


-1.68894e-02 


1.306706-02 


-2.559686-03 


df 


2.98626e+00 


2.57075e+01 


-3.24294e+00 


5.579426+00 


-3.346656+00 


df 


2.70854e+00 


2.82522e+01 


1.82688e+01 


-1.703786+01 


2.777756+00 


df 


7.22072e-01 


2.16565e+01 


4.64775e+00 


-5.30941C+00 


8.56825e-01 


df 


-5.90865e-02 


4.25941e+00 


-2.83516e+00 


1.866896+00 


-3.960366-01 








DBW2 






c« 


4.47375e-02 


5.75118e-04 


-3.46013e-04 


8.997876-05 


9.519946-06 




1.91634e-03 


-2.99489e-02 


2.04894e-02 


-3.254906-02 


1.009406-02 


(«) 
"2 


4.13973C-02 


9.03824e-02 


-2.27829C-01 


2.263706-01 


-8.34536C-02 


(«) 
"3 


6.61204e-02 


2.46242e-01 


2.16915e-01 


-8.469066-02 


-9.705256-03 




-3.73368e-03 


6.80095e-02 


-6.85793e-02 


5.411156-02 


-1.249666-02 


df 


1.32735e+01 


5.46203e+01 


-2.94128e+01 


7.785186+01 


-3.162316+01 


df 


2.63415e+01 


8.32744e+01 


6.80241e+01 


-5.891696+00 


-1.194116+01 


df 


1.15204e+01 


6.95572e+01 


2.59625e+01 


-5.711566+00 


-3.67863e+00 


df 


-7.17896e-01 


1.26939e+01 


-1.26694e+01 


9.674106+00 


-2.265546+00 



TABLE I: Numerical results of the fit parameters for A. 
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i 









1 




2 




3 




4 












plaquette 












4.52373e- 


-02 


7.95780e- 


-06 


5.08306e- 


-04 


-4.22552e- 


-04 


6.990056- 


-05 


(») 
n\ ' 


2.84911e- 


-03 


7.82834e- 


-04 


-3.98343e- 


-03 


7. 842536- 


-04 


1.646356- 


-03 


(i) 

n\' 


1.04294e- 


-02 


-3.50230e- 


-03 


2.54349e- 


■02 


-8.17546e- 


-03 


2.133006- 


-03 




-4.07134e- 


-03 


1.79561e- 


-02 


-2.42208c- 


02 


1.70135e- 


-02 


-4.40267c- 


-03 


nf 


6.34372e- 


-04 


-4.76941e- 


-03 


1.51021e- 


-02 


-6.08709e- 


-03 


-6.67261e- 


-05 


df 


1.06023e+00 


1.31409e+00 


1.05189e+00 


7.46079e- 


-01 


2.098686- 


-01 


df 


1.06865e+00 


8.17638e- 


-01 


7.32570e- 


-01 


4.22780e- 


-01 


-5.17816e- 


-01 


df 


-4.28339e- 


-01 


9.00811c- 


-01 


8.22650c- 


-01 


3.40122e- 


-01 


-4.14352c- 


-01 


df 


6.74709e- 


-02 


-4.02057e- 


-01 


1.18565e+00 


-7.43252e- 


-01 


1.207466- 


-01 












Iwasaki 














4.49640e- 


-02 


1.43870e- 


-04 


1.38586e- 


-04 


-8.90746e- 


-05 


-4.85512e- 


-06 


(») 
n\ ' 


2.55443e- 


-03 


-1.78947e- 


-03 


-1.16834e- 


-03 


7.84991e- 


-06 


6.304846- 


-04 


(j) 


6.56198e- 


-03 


4.01176c- 


-03 


9.98744c- 


-03 


-3.19040c- 


-04 


-9.73106c- 


-04 


(«) 
n\' 


-2.32566e- 


-03 


8.53741e- 


-03 


-9.25619e- 


■03 


6. 135656- 


-03 


-2.25368e- 


-03 


nf 


2.48873e- 


-04 


-1.39051e- 


-03 


5.97769e- 


■03 


-2.99920e- 


-03 


2.505806- 


-04 


df 


6.08980e- 


-01 


1.33407e+00 


9.80322e- 


■01 


5.334466- 


-01 


-1.03361e- 


-01 


df 


8.28020e- 


-01 


8.06319e- 


-01 


6.73518e- 


-01 


3.080316- 


-01 


-6.20501e- 


-01 


df 


-3.09041e- 


-01 


6.39750c- 


-01 


2.78380e- 


-01 


-1.60949c- 


-01 


-2.75471c- 


-02 


df 


3.30870e- 


-02 


-1.23565e- 


-01 


4.61016e- 


■01 


-2.91189e- 


-01 


4.433516- 


-02 












DBW2 












c^) 


4.49949e- 


-02 


-5.33912e- 


-05 


9.92059e- 


■04 


-9.12758e- 


-04 


2.316026- 


-04 


n\ ' 


1.25755e- 


-05 


-2.32660e- 


-03 


-9.17352e- 


■03 


1.000976- 


-02 


-3.13132e- 


-03 


(i) 


1.17640c- 


-02 


7.59642c- 


-03 


8.23216c- 


-03 


-9.46580c- 


-03 


4.57219c- 


-03 


n\' 


-2.35282e- 


-03 


4.99363e- 


-03 


-3.60508e- 


■03 


9.554216- 


-03 


-5.03662e- 


-03 


nf 


7.79286e- 


-05 


5.07874e- 


-04 


7.65744e- 


-04 


-1.54716e- 


-03 


5.609276- 


-04 


df 


8.38142e- 


-01 


9.34378e- 


-01 


6.79401e- 


-01 


4.382166- 


-01 


1.826176- 


-01 


df 


1.79736e+00 


1.35771e+00 


9.74925e- 


-01 


3.911386- 


-01 


-7.30343e- 


-01 


df 


-4.31666e- 


-01 


3.18664c- 


-01 


1.50605c- 


-01 


1.16987e- 


-01 


-1.26641c- 


-01 


df 


2.39747e- 


-02 


8.19427e- 


-02 


6.10548e- 


-04 


-8.72171e- 


-02 


3.347826- 


-02 



TABLE II: Numerical results of the fit parameters for A. 
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i 





1 


2 


3 


4 








plaquette 








9.84230e-03 


-3.51914e-03 


6.52579e-03 


-2.364386-03 


3.701716-04 


(») 
n\ ' 


2.78548e-02 


-4.71051e-01 


7.97658e-01 


-1.588306+00 


5.932286-01 


(i) 

n\' 


5.84173e-01 


1.40053e-01 


7.76400e-01 


-2.644566+00 


9.174756-01 


n\> 


-1.26248e+00 


-1.11733e+00 


-3.079806+00 


7.89903C-01 


6.066796-01 


nf 


1.52550e-01 


-1.26536e+00 


1.14001e+00 


-9.708846-01 


2.542406-01 


df 


2.24515e+01 


1.27087e+02 


1.13085e+02 


-6.903816+01 


-6.791946+00 


df 


5.48959e+01 


1.19734e+02 


2.07704e+02 


-1.603316+02 


2.149156+01 


df 


1.05293e+02 


1.66311C+02 


2.45242e+02 


-2.242236+02 


2.135626+01 


df 


-1.42301e+01 


1.24790e+02 


-9.07259e+01 


5.730376+01 


-1.527926+01 








Iwasaki 






c(^) 


7.10943e-03 


-8.49976e-03 


1.02419e-02 


-5.117176-03 


1.091726-03 




9.18767e-03 


-7.17314e-02 


2.20257e-01 


-2.187206-01 


5.116346-02 


n\' 


4.32270e-02 


7.01214e-02 


-1.338666-01 


5.226296-02 


-3.098626-02 


n\' 


-1.61351e-01 


-1.04229e-01 


-8.888776-03 


-2.639906-01 


1.153326-01 


nf 


5.91548e-03 


-3.94141e-02 


5.92655e-02 


-9.094456-02 


2.175756-02 


df 


3.12332e+00 


1.55819e+01 


1.37894e+01 


1.118196+01 


-1.025386+01 


df 


8.14387e+00 


4.36052e+00 


1.525286+01 


7.734706+00 


-3.223516+00 


df 


1.83424C+01 


2.00024e+01 


1.545356+01 


9.908106+00 


-1.026026+01 


df 


-5.89261e-01 


5.98743e+00 


-7.246676+00 


1.553186+01 


-5.126796+00 








DBW2 






c^) 


3.24314e-03 


2.52054e-03 


4.444316-03 


-3.084516-03 


6.313746-04 


n\ ' 


1.05779e-01 


-1.01997e+00 


1.350986+00 


-2.362836+00 


9.061066-01 


n\' 


1.34493C-01 


7.58544C-01 


-1.696566+00 


6.160426-01 


-1.907206-01 


n\' 


-7.86597e-01 


-2.38384e+00 


4.729236-03 


-2.825196+00 


1.725446+00 


nf 


6.41759e-02 


-2.48777e-01 


-6.187216-02 


l.lOMle-Ol 


-5.586596-02 


df 


4.67326e+01 


6.66294e+01 


9.217886+00 


1.056576+02 


-6.267516+01 


df 


-1.58592e+01 


7.51684e+01 


3.202306+01 


-4.757316+01 


2.294016+01 


df 


1.66161C+02 


1.31914C+02 


1.135636+02 


4.161936+01 


-7.633906+01 


df 


-1.23813e+01 


8.96851e+01 


-8.013926+01 


3.673806+01 


-5.422576+00 



TABLE III: Numerical results of the fit parameters for Cy^ . 



22 



i 





1 


2 


3 


4 








plaquette 






c(») 


4.72138e-03 


1.18712e-02 


-3.70767e-02 


3.644176-02 


-1.057576-02 


(») 
n\' 


2.88129e-01 


-1.02535e+00 


2.46673e+00 


-2.18206e+00 


5.911466-01 


n\' 


1.94450e-01 


6.18542e-02 


5.07036e-02 


-2.53440e-01 


6.961766-02 


(i) 


1.12189e-01 


-7.65970e-01 


2.087636+00 


-1.879716+00 


5.163206-01 


(«) 

nl' 


6.00100e-02 


5.83904e-02 


5.29609e-03 


-7.12338e-02 


3.055376-02 


df 


2.63902e+01 


1.69601e+01 


5.43255e+00 


-1.02071e+01 


4.196826+00 


df 


1.94234e+01 


2.93848e+01 


1.60899e+00 


-2.29237e+01 


1.077806+01 


df 


1.26822e+01 


3.48188C+01 


5.02924e+00 


-3.229976+01 


1.47403C+01 


df 


2.76062e+00 


1.15866e+01 


-1.75853e+01 


1.099956+01 


-2.314016+00 








Iwasaki 






c« 


1.53152e-03 


1.03025e-02 


-3.42414e-02 


3.466776-02 


-1.003686-02 




3.09889e-01 


-1.04777e+00 


2.55683e+00 


-2.288076+00 


6.099056-01 


(i) 


-6.99951e-02 


-4.40763e-01 


1.04528C+00 


-1.079996+00 


2.662736-01 


(«) 
n\' 


4.81714e-02 


-6.08612e-01 


1.06819e+00 


-7.596636-01 


1.486896-01 


(*) 


7.12031e-03 


-4.51084e-02 


1.18638e-01 


-1.148916-01 


3.342776-02 




3.28468e+01 


2.10820e+01 


-1.72327e+00 


-4.859046+00 


2.967236+00 


d« 


2.13491e+01 


3.82549e+01 


1.27385e+01 


-3.262906+01 


1.424846+01 


df 


7.91264e+00 


3.20460e+01 


1.796116+00 


-3.172876+01 


1.562336+01 


df 


1.42775e+00 


5.13906e+00 


-9.29549e+00 


7.444616+00 


-2.044736+00 








DBW2 






c« 


3.59477e-03 


-2.50302e-03 


-7.48782e-03 


1.208536-02 


-4.332866-03 




-4.65181e-03 


-1.36855e-01 


3.58241e-01 


-3.410296-01 


1.014826-01 


(*) 

"2 


-1.55086C-01 


-1.09153C-01 


2.479966-02 


-3.567616-02 


-1.714466-02 


(») 
"3 


6.09454e-03 


-1.03449e-01 


1.73653e-01 


-1.565216-01 


5.348036-02 


(*) 

n\' 


-6.32341e-03 


-1.77005e-01 


3.75656e-02 


1.072676-01 


-5.561196-02 


df 


5.18966e+00 


2.00129e+00 


3.24520e-01 


1.372216+00 


-1.015916-01 


df 


6.12021e+00 


7.41829e+00 


1.22181e+00 


-4.718976+00 


3.050916+00 


df 


1.89231e-01 


4.22273e+00 


-6.553296-02 


-2.043066+00 


5.660556-01 


df 


3.67961e-01 


6.56838e+00 


9.83657e-01 


-7.059326+00 


3.114986+00 



TABLE IV: Numerical results of the fit parameters for c^ 
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i 







1 




2 




3 




4 










plaquette 












-3.77901e-03 


-5.20628e- 


-03 


6.92389e- 


-03 


-3.33440e- 


-03 


5.83119e- 


-04 


(») 
n\ ' 


-5.34436e-05 


1.86231e- 


-02 


-3.55829e- 


-02 


2.01514e- 


-02 


-5.51715e- 


-03 




4.10664e-03 


-1.02295e- 


-02 


6.70875e- 


■02 


-5.09115e- 


-02 


1.34993e- 


-02 


ny 


-7.67466e-03 


1.17796e- 


-02 


-3.73396e- 


02 


2.86610c- 


-02 


-5.86918c- 


-03 


nf 


7.17646e-03 


4.63334e- 


-03 


3.45373e- 


-03 


-2.78005e- 


-03 


-6.70178e- 


-04 


df 


6.73021e-01 


1.34845e+00 


1.32647e+00 


9.83849e- 


-01 


-6.91234e- 


-01 


df 


8.11070e-01 


1.79291e+00 


1.05732e+00 


4.24515e- 


-01 


1.04878e- 


-01 


df 


-2.10312C-01 


-3.58696c- 


-01 


-2.70077e- 


01 


-8.84022c- 


-02 


1.93194c- 


-01 


df 


1.54957e+00 


2.32307e- 


-01 


5.35889e- 


-01 


5.32861e- 


-01 


-5.67951e- 


-01 










Iwasaki 














-8.92573e-04 


-4.98763e- 


-03 


5.83872e- 


-03 


-2.80997e- 


-03 


4.99720e- 


-04 


n\ ' 


8.18233e-03 


1.46155e- 


-02 


-2.51372e- 


-02 


1.74577e- 


-02 


-5.53396e- 


-03 


ny 


-1.08859e-02 


-4.35439e- 


-03 


4.12491c- 


02 


-3.05207c- 


-02 


1.01006e- 


-02 


ny 


7.18333e-03 


1.37046e- 


-02 


-2.22405e- 


■02 


1.97996e- 


-02 


-6.71756e- 


-03 


nf 


3.10051e-03 


3.88076e- 


-03 


-4.07139e- 


■03 


1.22677e- 


-03 


-3.58180e- 


-04 


df 


2.58979e-01 


4.42756e- 


-01 


6.74009e- 


-01 


6.44566e- 


-01 


-4.79794e- 


-01 


df 


1.05613e+00 


1.31787e+00 


9.98757e- 


-01 


6.30569e- 


-01 


1.67564e- 


-02 


df 


4.35593e-01 


6.79640e- 


-01 


5.08309e- 


-01 


2.20645e- 


-01 


-2.49549c- 


-01 


df 


1.68734e+00 


-4.31052e- 


-01 


1.06194e- 


-01 


2.17586e- 


-01 


-2.29057e- 


-01 










DBW2 












c^) 


1.65353e-03 


-5.20018e- 


-03 


5.58906e- 


■03 


-2.68645e- 


-03 


4.85235e- 


-04 


n\ ' 


3.05240e-02 


1.73659e- 


-02 


-2.07260e- 


■02 


1.51523e- 


-02 


-5.59690e- 


-03 


(i) 

ny 


-3.78827e-02 


1.51961e- 


-02 


1.31476c- 


-02 


-9.24480c- 


-03 


4.14817c- 


-03 


ny 


3.28025e-02 


6.55780e- 


-03 


-1.23069e- 


■02 


1.14891e- 


-02 


-6.25979e- 


-03 


nf 


-1.50078e-03 


1.08107e- 


-02 


-9.76440e- 


■03 


2.91846e- 


-03 


-2.97660e- 


-04 


df 


1.69736e+00 


6.54862e- 


-01 


6.24044e- 


-01 


4.66231e- 


-01 


-4.83184e- 


-01 


df 


1.75714e+00 


6.50739e- 


-01 


4.07680e- 


-01 


2.75134e- 


-01 


-7.55484e- 


-02 


df 


8.31723e-02 


1.05494e+00 


8.70566c- 


01 


3.58793e- 


-01 


-4.81543e- 


-01 


df 


2.94855e+00 


-1.89555e- 


-01 


2.11226e- 


-02 


-1.18785e- 


-01 


-1.54914e- 


-01 



TABLE V: Numerical results of the fit parameters for c^ 
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i 





1 


2 


3 


4 








plaquette 






c(») 


3.50525e-03 


-4.68207e-02 


8.77916e-02 


-6.22507e-02 


1.481696-02 


(») 
n\' 


-3.76844e-01 


2.12957e+00 


-4.14093e+00 


3.003546+00 


-7.257576- 01 


n\' 


-2.06844e-01 


7.70017e-01 


-9.00837e-01 


4.254476-01 


-7.252606-02 


n\' 


-2.78087e-01 


2.41752e+00 


-4.31448e+00 


2.949446+00 


-6.346836-01 


n\' 


-1.82981e-01 


2.78762e-01 


-4.96115e-01 


3.481066-01 


-9.336426-02 


df 


3.68004e+01 


7.08098e+00 


-3.92992e-02 


-5.208576+00 


1.473756+00 


df 


1.30394e+01 


2.67925e+01 


-6.61332e+00 


-2.847796+01 


1.771526+01 


df 


3.07591e+01 


2.40240e+01 


-1.04165C+00 


-1.790106+01 


1.016106+01 


df 


5.72652e+00 


1.12568e+01 


-1.88964e+01 


1.174696+01 


-1.684786+00 








Iwasaki 






c« 


1.53487e-03 


-3.41468e-02 


6.82220e-02 


-5.090286-02 


1.277286-02 




-8.27139e-02 


1.34471e+00 


-2.86275e+00 


2.207716+00 


-5.682656-01 


(«) 


1.03514e-01 


1.49986e+00 


-2.05986e+00 


1.135426+00 


-2.110336-01 


n\' 


1.31475e-02 


2.05314e+00 


-3.25439e+00 


2.047006+00 


-4.107816-01 




-4.37623e-02 


2.63972e-01 


-4.70760e-01 


3.266376-01 


-7.590316-02 




3.30552e+01 


4.09568e+00 


-5.98710e+00 


6.093426+00 


-2.340396+00 




1.77638e+01 


4.26837e+01 


-6.91271e+00 


-3.911596+01 


2.027906+01 


df 


3.03378e+01 


2.33940e+01 


-4.557826+00 


-1.45824e+01 


7.893106+00 


df 


3.68914e+00 


8.68953e+00 


-1.06404e+01 


3.382076+00 


6.353096-01 








DBW2 






c« 


2.37474e-03 


-1.26844e-02 


2.50575e-02 


-2.076706-02 


5.766656-03 




6.48432e-02 


1.36343e-01 


-3.45284e-01 


3.370776-01 


-1.043086-01 


(*) 
n\' 


5.92327C-02 


1.00062e-01 


-5.62452e-02 


2.127276-02 


1.201386-02 


(») 


1.04664e-01 


4.89947e-01 


-4.13500e-01 


1.426956-01 


2.248286-03 


n\' 


-1.51902e-03 


1.41416e-01 


-6.39386e-02 


-1.099496-01 


6.863656-02 


df 


4.96718e+00 


-1.30295e+00 


-1.39159e+00 


6.265646+00 


-2.707256+00 


df 


2.70372e+00 


8.98447e+00 


-1.36186e+00 


-8.099356+00 


4.706016+00 




6.12393e+00 


8.34213e+00 


-5.83018e-01 


-3.189816+00 


1.538686+00 




1.28638e-01 


6.27763e+00 


-3.24311e+00 


-4.434416+00 


2.887246+00 



TABLE VI: Numerical results of the fit parameters for ci 
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i 





1 


2 


3 


4 








plaquette 






c(») 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


(») 
n\' 


2.09397e-03 


-5.33186e-04 


-4.65192e-03 


5.42610e-03 


-2.24704e-03 


n\' 


-1.06471e-02 


2.27643e-02 


4.01977e-02 


-6.10584e-02 


2.45592e-02 


n\' 


1.38706e-02 


5.94564e-02 


1.64618e-02 


4.22552e-02 


-3.488156-02 


n\' 


-1.08445e-03 


-1.93452e-03 


7.00296e-03 


-8.09199e-03 


2.98526e-03 


df 


1.22751e+00 


2.94194e+01 


1.07318e+01 


-5.53447e+00 


-1.88793e+00 


df 


-9.43352e+00 


1.44250e+01 


7.73519e-01 


-6.89522e+00 


3.22929e+00 


df 


1.12565e+01 


1.83971e+01 


1.66376C+01 


-1.04853C+00 


-7.23085C+00 


df 


-1.75364e-01 


2.06473e+00 


8.42623e+00 


-7.65891e+00 


1.65849e+00 








Iwasaki 






c« 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


(») 

ny 


8.71136e-03 


1.34446e-02 


-6.72308e-02 


6.01238e-02 


-1.63627e-02 


(«) 


-4.80524e-02 


9.21777C-02 


5.29502C-01 


-6.090206-01 


1.71269e-01 


n\' 


5.71428e-02 


5.68763e-01 


-3.91779e-01 


2.55005e-01 


-1.056006-01 


(») 


-6.84520e-03 


-7.56285e-03 


1.02776e-02 


-l.llOOOe-02 


5.96372e-03 


d« 


2.36201e+01 


1.85591e+02 


1.27770e+02 


-2.20287e+02 


5.36880e+01 


d« 


-7.63923e+01 


8.29550e+01 


-7.43883e+01 


5.61588e+01 


-9.224756+00 


df 


6.70337e+01 


2.16961C+02 


-4.95764e+01 


-4.63254e+01 


1.00042e+00 


df 


-4.14600e+00 


4.16800e+01 


-4.09622e+01 


2.30321e+01 


-5.120056+00 








DBW2 






c« 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 




1.24968e-02 


9.44758e-03 


-3.42722e-02 


3.14005e-02 


-9.984436-03 


(*) 
n\' 


-3.82921e-02 


1.39231C-01 


4.08137e-01 


-5.22501e-01 


1.54012C-01 


(i) 


4.77557e-02 


4.30929e-01 


-3.25076e-01 


2.27980e-01 


-9.176446-02 


(») 
n\' 


-5.77412e-03 


-1.55592e-02 


1.13224e-02 


-8.06948e-03 


4.26375e-03 


df 


3.83004e+01 


1.21779e+02 


1.19599e+02 


-1.84372e+02 


4.54188e+01 


df 


-9.12407e+01 


5.90943e+01 


-7.77565e+01 


5.45309e+01 


-4.401116+00 


4) 


6.95720e+01 


2.16484e+02 


-4.76922e+01 


-4.190546+01 


-8.566106-01 




-3.56676e+00 


2.24451e+01 


-2.61850e+01 


1.87052e+01 


-4.820866+00 



TABLE VII: Numerical results of the fit parameters for 
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i 





1 


2 


3 


4 








plaquette 






c(») 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


(») 
n\' 


-4.92854e-02 


7.76273e-02 


-1.09117e-01 


5.41335e-02 


-6.65490e-03 




1.62369e-02 


-2.05156e-01 


4.87357e-02 


-1.00512e-02 


2.19379e-02 


n\' 


-1.62807e-02 


-6.000906-02 


9.26995e-03 


-9.52634e-02 


5.14551e-02 


(i) 

n\' 


-1.34504e-02 


-1.09913e-01 


-9.93819e-03 


3.72108e-02 


-1.88461e-03 


df 


1.07320e+00 


8.61247e+00 


-1.70831e+00 


2.57570e+00 


-2.35960e+00 


df 


-6.07917e-01 


7.62319e+00 


4.48235e+00 


1.52722e+00 


-2.26997e+00 


df 


8.80495C-01 


8.95034e+00 


1.54057e+00 


-8.81352e-02 


-1.36898C+00 


df 


2.21861e-01 


2.61262e+00 


5.61100e-01 


-1.25562e+00 


1.61500e-01 








Iwasaki 






c« 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


(») 
n\' 


-1.89216e-02 


7.61233e-02 


-9.73709e-02 


5.06130e-02 


-8.17803e-03 




-7.31825e-03 


-1.74543e-01 


2.94655e-01 


-1.91881e-01 


5.14398e-02 


(«) 
n\' 


-7.05824e-02 


2.71060e-01 


-1.02481e-01 


1.73657e-01 


-6.37448e-02 


n\' 


-9.00773e-03 


-1.67444e-01 


1.36274e-01 


-1.33701e-01 


5.65792e-02 


df 


-1.97613e+00 


3.89137e+01 


1.57804e+00 


2.54447e+00 


-6.67567e+00 


df 


4.57630e+00 


3.62223e+01 


2.26543e+01 


9.04292e+00 


-1.38239e+01 


df 


3.86975e+00 


3.56956e+01 


3.56635e+01 


1.22621e+01 


-1.85228e+01 


df 


8.34828e-02 


1.73013e+01 


7.67351e+00 


-3.68814e+00 


-1.92297e+00 








DBW2 






c« 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


O.OOOOOe+00 


(») 
n\' 


3.94399e-02 


6.12762e-02 


-4.42999e-02 


2.91012e-03 


5.62763e-03 




-5.50451C-04 


1.53894e-02 


7.33107e-02 


-2.14443C-02 


1.26660e-02 


ny 


1.44008e-02 


-5.31159e-03 


5.91542e-02 


-7.08047e-02 


3.46751e-02 


(*) 
n\' 


2.32681e-02 


8.14410e-02 


-3.26852e-02 


6.02456e-02 


-3.77997e-02 


df 


1.45225e+00 


1.06854e+00 


8.00493e-01 


4.90831e-01 


-1.84695e-01 


df 


4.10400e-01 


6.60514e-01 


6.29942e-01 


5.51064e-01 


2.32655e-01 


4) 


1.24019e+00 


1.09508C+00 


7.99430e-01 


4.32776e-01 


-4.18890C-01 


df 


8.12694e-01 


1.16732e+00 


7.68715e-01 


2.41916e-01 


-6.06008e-01 



TABLE VIII: Numerical results of the fit parameters for ch 
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FIG. 1: One-loop diagrams for the vertex functions, q denotes the outgoing momentum and p denotes 
the incoming momentum. 
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0.044 
0.06 
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